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Morris Perlmutter 

Lewis Research Center 
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•Abstract 

Couette flow and heat transfer through a hard-sphere 
rarefied gas enclosed between parallel walls are analyzed 
by a Monte Carlo procedure. In one case, both walls are 
stationary and at different temperatures. In the other 
case, the upper wall is moving, and both walls are at the 
same temperature. The molecules are assumed to have hard 
sphere collision and a wall accommodation factor of unity. 
The wall temperature ratio and wall velocity are taMen large 
so that linearized solutions are not applicable. Sample 
molecule histories are followed across the channel, which 
has been divided into zones. The target molecules in each 
zone are assumed to have a different two-sided half- 
Maxwellian velocity distribution. By scoring the properties 
of the sample molecules as they pass scoring positions, the 
macroscopic quantities of interest, such as temperature, 
density, shear, and heat transfer, are obtained. Also, the 
average thermal velocity and number density in each zone 
are found and compared with the values assumed in the target 
molecule distribution. The problem is iterated until these 
values agree. 


Symbols 

C = average thermal velocity (2RgT)-*-A 
D - channel width 

f = probability' distribution function 
Kn = Knudsen number Aj^/D, = m /-v/2 p^Dva^ 
M = dimensionless velocity u/c^-' q : 
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RAREFIED GAS DYNAMICS 
= mass of molecule 

= number of sample molecules emitted from surface 0 
in Monte Carlo run, proportional to flux of mole- 
cules leaving surface 0 
= pressure; Pii/3 
= shear stress (-p(V^V-)) 

= zone or scoring -position number 
= last scoring-position number 
= property of sample molecule 

- averaged quantity Jo,f d-% 

= averaged quantity yQf + d 3 V 

= heat transfer 
= total energy transfer 
= random number between 0 and 
= gas constant 

= mutual collision cross section 

= number of sample molecules through scoring position 
in positive, negative xg direction 
= absolute temperature 
= dimensionless temperatures T/T w q 
= mean velocity 
= coordinate 
= molecular velocity 
= thermal velocity; v - u 
= velocity after collision 
= defined by equations (B5) and (B6) 

= defined by equations (B5) and (B6) 

= viscosity coefficient 
= thermal conductivity 

= collision rate of sample molecule with target mole- 
cules 

= path length to collision 

= mean free path length through a Maxwellian gas with 
density p and thermal velocity C 
= mean free path length 
= dimensionless velocity, V'/C 
= defined by equation (All) 

= density 

= diameter of hard sphere molecule 
= angle between sample molecule and target molecule 
velocities 


Subscripts : 

fl = averaged across channel 
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RAREFIED GAS DYNAMICS 
= Maxwellian 

= refers to case with Doth walls stationary 
= increment number 
= last increment 
= relative velocity 

= radial coordinate in cylindrical system 
= sample molecule 
= target molecule 
w,0;w, 1 = at wall 0,1 
0,1 = evaluated next to 0,1 

1,2,3 = coordinate directions 

+,- = positive or negative Xg direction 

Introduction 

The present -analysis treats a rarefied gas enclosed be- 
tween two walls. Two cases are treated: in one, both walls 

are stationary and at different temperatures, while in the 
other, both walls are at the same temperature but one wall 
is moving. The molecule are assumed to have hard-sphere 
collisions and a wall accommodation of 1. The problem is 
analyzed .by .-a Monte Carlo technique, that is, a model 
sampling procedure. Sample molecule histories are gener- 
ated by making -choices at points of decision from the appro- 
priate probability distributions. These sample molecules 
are followed through the model and, by averaging certain 
molecular properties at various positions, the macroscopic 
quantities of interest can be obtained. This method reduces 
the complexity of the analysis since the relations for a 
single sample molecule history are relatively simple. Thus, 
simplifying assumptions made in the more usual analytical 
procedures are also avoided. However, the Monte Carlo 
method requires extensive calculations on a high-speed com- 
puter. for the many sample histories needed. 

Gross and.Ziering (ref. l) treated the couette flow 
problem for the' case of low relative plate velocities with 
both plates at the same temperature. They also treated the 
case where both walls were stationary but at different tem- 
peratures (ref. 2). Because small differences in relative 
wall velocities and wall temperature were assumed in these 
analyses, the linearized form of the Boltzmann equation 
was used. The local distribution function was approxi- 
mated -by half -range -polynomials ''in velocity space.'-The'space- 
dependent' coefficients in this distribution were found by 
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RAREFIED GAS DYNAMICS 

forming half -range moment equations. 

Lui and Lees (ref. 3) have analyzed the ease of rar- 
efied couette flow with heat transfer assuming Maxwellian 
molecule collision. They assumed a local "two -sided half- 
Maxwellian velocity distribution. This is substituted into 
the moment equations , which are then solved for 'the parame- 
ters of the distribution. From this, the macroscopic floT 
parameters can "be found. 

This same approach was used by Lavin and Haviland 
(ref. 4) for the case of stationary walls at different te' 
peratures assuming hard sphere molecule collisions. 

The Monte Carlo method has been used previously in 
rarefied gas problems 'by Perlmutter (ref. 5) to treat a 
collisionless plasma flowing -through a channel with an im- 
posed magnetic field. Haviland (ref. 6) used Monte Carlo 
to treat a rarefied gas with hard-sphere collision con- 
tained between stationary parallel plates at different wa! 
temperatures. In that analysis, sample molecule histories 
are followed through collision with target molecules as 
'they -travel through -the model. From the histories of the 
sample molecules, tables of the distribution of the veloc- 
ity components can be obtained in the various zones into 
which the channel has been divided. These tables can then 
be used for the target molecule velocity distribution in 
■the next iteration of the sample molecule histories. Aftex 
the converged distribution of velocities are obtained, the 
macroscopic quantities of interest, such as temperature, 
density, and heat transfer, can be obtained by finding the 
moments of the distribution in each zone. 

The present analysis also follows a sample molecule 
traveling through the model. However , in the present "-case, 
the target molecule velocities in each zone -are assumed to 
have a two-sided half -Maxwellian distribution of the form - 
'assumed by Lui and Lees in reference 3. By scoring proper- 
ties of the sample molecules at scoring position across the 
channel, the parameters of the distribution as well as the 
macroscopic flow properties such as- temperature, velocity 
profile, heat transfer, and shear stress can be found. This 
method was used by Perlmutter (refs. 7 and 8) for heat 
transfer between parallel plated with the walls stationary 
but at different temperatures. In that solution, it was 
assumed that the local distribution of target molecules was. 
Maxwellian instead of two-sided half -Maxwellian. 
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Sample History 


The model analyzed is shown in figure 1. The wall at 
Xg of 0 is stationary and at temperature T w q. The 
other wall at Xg of 1 is moving at velocity’ u^. ^ and 
is at temperature T w The distance across the channel is 
divided into zones, in each zone, the target molecule ther- 
mal velocity distribution is assumed to he a two-sided half- 
Maxwellian distribution given by 



( 1 ) 


where f+ refers to the molecules moving in the positive 
Vg direction and f_ refers to molecules moving in the 
negative V 2 direction. 

It is assumed that the molecules reflected are per- 
fectly accommodated, that is, in a Maxwellian distribution 
at the wall based on the wall temperature. The velocity of 
the sample molecule leaving the wall as derived in refs. 6 
and 8 is given by 

v 2 = (-eg m % 2 ) 1/2 (2a) 


Vn 


= (-eg 


ln % r ) 


l/2 


cos 2itR q + up w 


(2b) 


V 3 = (-eg In R Vr ) 1 / 2 sin 2*R§ 


(2c) 


where R is a random number between 0 and 1 chosen for 
each' sample molecule leaving the surface. 

Path length to collision . - After the sample molecule 
leaves the wall, its path length to collision must be calcu- 
lated to find if it has a collision before passing through 
the first zone. -The probability that a sample molecule will 
collide in the incremental path length A to A + dA is 
given (refs. 8 and 9) as 

exp(-A/Aj 

AT" (3) 

where A s is the mean free path to collision of the -sample 
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molecule moving -at velocity V s in that -zone. A path 
length to collision for the sample molecule can he chosen 
from this distribution- hy 'the same procedure of picking a 
random number used earlier -for the sample molecules leaving 
the surface: 

A = -A s In R A (4) 

To pick from this distribution, the mean free path A s for 
the sample molecule in the zone must be known. If the un- 
directed velocities of the target molecules are in a Max- 
wellian distribution with, an average thermal velocity C 
and density p, ' the mean free path for the sample molecule 
Ajj> as shown in appendix A, is given by 



where Kn is the usual definition of ‘Khudsen number for a 
gas in a Maxwellian distribution with hard-sphere molecule 
collisions. Since the target molecule velocity distribu- 
tion has been assumed to be a two-sided half -Maxwellian 
distribution (eq. (l) ) , the mean free path is obtained by 
assuming the distribution of target molecules to consist of 
two full Maxwellian distributions in each increment 


f t = %f + f M- 


( 6 ) 


where f^ + 'and %_ are Maxwellian, distributions with deri- 
■siti,es 2 pi and 2p_ and thermal velocities of C + and 
C,_, respectively. However, the molecules of fj^ + with 
Vg < 0 and the- molecules of f M- with Vg > 0 are con- 
sidered nonexistent. This permits the calculation of the 
mean free path using the Maxwellian distribution relations. 
Thus, the mean free path for a -sample molecule moving 
through two Maxwellian gases, as shown in appendix A, is 
given by 

-A s A M+ A m _ 

where Ajyj-f. or Am_ is evaluated, as in equation (5). 

After a path length to collision is picked for the 
sample molecule from equation (4), it is compared with the 
-distance the sample molecule must travel to cross the zone. 
If A is greater than this distance, the sample molecule is 
started at the beginning of ‘the next zone with its -velocity 
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components unchanged and the procedure is as "before. If it 
is not, there is a collision in the zone- and new velocity 
components -must "be calculated for the sample molecule after 
collision. 

Sample Molecule Velocity After Collision 

To find the new sample molecule velocity -components 
after collision, the velocity components of the target 
molecule collision partner -must first "be found. Since there 
are two ■ distributions of target molecules locally, fjj + and 
%-> the distribution the target molecule collision partner 
is to be picked from must first be decided. The fraction of 
all the collision partners that come from the fj,j+ distri- 
bution is given by A s /A$/[+* Then, for -a picked random num- 
ber, ^.f R < A s /Am+.> the collision partner velocity compo- 
nents are picked from the f M+ distribution. Otherwise, 
they are picked from ’the %_ distribution. The equations 
for obtaining the velocity components of the target molecule 
collision partner are given in appendix 

If Vg-fc picked from the f^ + distribution is negative 
or, if picked from fjji, if V 2 t is positive, the collision 
partner is considered nonexistent. This is because half the 
molecules in the f M+ and %_ distribution are nonexist- 
ent, as discussed before. In this case, the sample molecule 
continues its history from the point of ‘collision with, its 
original velocity. 

After collision of the sample molecule with the chosen 
collision partner, the new components of velocity V* of 
the sample molecule are given (ref. 8) by 

’is - 1 (As + At) + V 1 - b2 ) l/2 H 1 < 8a ) 

As = I (As + At) + At 1 - 21,2 ) > ( 8b ) 

As = | (As + At) + Ad - b2 ) i/2 3 J < 8c ) 

where Vr -is given by equation (A2). The values of H and 
E are, obtained by -picking two random numbers which are used 
in the following equations: 

- H = 2R h - 1, E = SR# - !, and b^ = H 2 + 'a (9) 

where b^ must be less than 1, or -a new set of random num- 
bers must be chosen to find H and S. 

The sample molecule history is then continued from the 
point of collision with its new velocity components. The 
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sample histories are completed -when the sample molecule re- 
turns to the surface from which it has "been emitted. 

Scoring -to Find Macroscopic Flow Properties 

The macroscopic fluid characteristics needed are den- 
sity, temperature, velocity, and heat transfer across the 
channel. These properties are obtained at scoring positions 
located at various distances across the -channel, as shown in 
figure 1. The average quantity of Q transported across 
the scoring -cross ■ section p, as shown in ref. 8, is 

p(v 2 Q> p =- (p + (V 2 Q) + ,) p + (p_<V 2 Q>_) p 



If Q is taken as l/Vg, 



Thus, by summing the inverse V 2 -velocity of the sample 
molecules as they -cross the- scoring position the local 
density -of' the., molecule 'with -a pbsitiye, negative compo- •' 
nent Of V 2 , p + , p_, can be obtained as well as the local density 
p. The average density -across the channel is then obtained 
by averaging the density of all the scoring cross sections. 

! If Q is taken as "Vq/V 2 , e'q. 10 becomes 



Dividing by the appropriate densities obtained' previously 
•g-ives the local mean velocity for the molecules with V 2 >0, 
u+, and V 2 < 0, u_, and also the local mean velocity, u. 

The undirected kinetic energy per unit mass of a mon- 
atomic gas is equal to the degrees of freedom times (l/2)RgT. 
This can' be written as 

(ik T - 1 - 1 <^> - f 4 

Then from equation (10) is obtained 


( 13 ) 
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The thermal velocity of the molecules with Vg > 0 
given .by- 



2 1 ,+ 
" 3 ri2 


is 

(15) 


and similarly for the molecules with V2 <0 



The shear stress across the channel is given by 
P 12 = p{VjVg) . This can he shown to he equal to 



The total energy being transferred across the channel 
is given by _ 

(V 2 ) (18) 

The transfer of the undirected kinetic energy across the 

channel is ,q = (l/2)p{V£(V’ 2 ) ), so that equation (18) be- 
comes 

9. = 9.|ji - u l P 12 ( 1®) 

By using these equations, the macroscopic quantities of 

interest are found. The values of p + , p_, 'C , C_ are then 
compared with the assumed values used in the xarget molecule 
distribution, and the problem is iterated until they are in 
agreement. Also obtained at the same time from the -above". ' 
equations -'.are the, local temperature,' density, ' shear,- and- • 
energy transfer. 


p +,0^w,0 
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Results and Conclusions 

Stationary walls with .a temperature ratio t^ q -°f 
0.25 was treated first. The local gas temperature divided 
by the temperature of the hot wall is shown in figure .2. 

The local density divided by the average density across the 
channel is given in figure 3. In figure 4, the heat flux 
transferred divided by the heat flux for the collisionless 
case, as given in appendix C, are shown. 

The present results can be seen to be in good agree- 
ment with the results of another Monte Carlo solution by 
Haviland (ref. 6). The linearized solution of Gross and 
Ziering (ref. 2) gives very good agreement with the present 
results except for the temperature profile. In the linear- 
ized solution of ref. 2, the centerline temperature is 
taken midway between the wall temperatures. The present 
results have their centerline temperature significantly be- 
low this value for the Knuds en number of 2; hence, the 
large disagreement in the results. 

The nonlinear results of Liu and Lees (ref. 3) of the 
temperature and density are of the correct magnitude but of 
■a somewhat different gradient. The nonlinear heat transfer 
results fall somewhat above the other solutions. 

Also shown in 'the figures are the slip continuum' re- 
sults." These results are obtained by using the fluid tem- 
perature in the continuum equations (appendix D) . By using 
■the slip wall temperatures that gave the best. fit both to 
the temperature and density results of the present solu- 
tion, the results shown were obtained. However, when the „ 
heat “transfer was calculated using these same slip wall tem- 
peratures and the continuum thermal conductivity, the heat 
transfer was found to be considerable above the other solu- 
tions. Since the slip -solutions with the continuum conduc- 
tivity are only expected to apply close to continuum condi- 
tions, this result is to be -expected. 

The heat transfer -at various positions across the chan- 
nel is shown in figure 6. Theoretically, this should be 
constant at each point across the channel (appendix D) , al- 
though the present method does not impose this condition in 
its method of solution. This condition, however, is satis- 
fied in the present results, as can be seen from this fig-, 
ure. 

To check the confidence limits of the present results, 
the solution was carreid out -20 -times with 2500 sample 
molecules each time for Kn =0.5. The temperature and 
density had a value of 5, the 95 percent confidence inter- 
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val, of close to 1 percent of their value , while the value 
of 8 for the heat transfer across the channel, was about 
2 percent of its value. The derivations and equation used 
in obtaining 8 are given in appendix E. 

The second case treated is when both walls are at the 
same temperature, but one wall is moving at a wall Mach num- 
ber M„. = u wl /C w ]_ of 4. The velocity profile is shown in 
figure 6. The linearized solution (ref. l) is in poor 
agreement with the present result, as would be expected for 
the nonlinear large 'wall Mach number case. 

The temperature profile across the channel is shown in 
figure 7. There is a substantial rise in temperature for 
smaller Knudsen numbers because, as the strongly directed 
stream of molecules leaves the moving wall, it -encounters 
the undirected stream. This gives a lower -value for the 
average directed motion of the fluid. This decrease in 
directed kinetic energy is now included in the thermal mo- 
tion which gives the higher temperature. 

The shear stress for the channel, shown in figure 8, 
has been nondimensionalized by the shear stress for the 
collisionless case (eq. (C5)). The results fall above the 
linearized solution of reference 1. The total energy 
transfer across the channel nondimensionalized by the total 
energy transfer -for the collisionless case (eq. (C6)) gives 
the same values as the dimensionless shear stress for the 
cases.. studied. This is also shown in figure 8, where shear 
stress and total heat transfer across the channel are plot- 
ted. Again, from the microscopic laws of conservation and 
the moment equations, it can be shown that the shear stress 
and total heat transfer at various positions across the 
channel should" be constant. Although this condition is not 
imposed by the present method of solution, the results in 
figure 10 indicate that the shear stress and total energy, 
transfer are nearly constant at different positions across 
the channel. 

In conclusion, the Monte Carlo method gives reason-, 
able results for the present problem without the necessity 
of extreme assumptions. The major drawback to the method 
is the increasing amount of computer time necessary to carry 
out the calculations as the Knudsen number became' smaller. 

Some typical running times for 10,000 sample molecules 
on an IBM 7094 were: Kn = 30, 8 minutes; Kn = 2; 10 min- 

utes; Kn = 0.5, 15 minutes.- These running times would make 
results difficult to^obtain for - small Knudsen numbers with 
the present -procedure. Generally, three iterations of- 
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50,000 samples were used. 

The present method is very flexible and can readily be 
used in other transport problems not readily solved by the 
more usual procedures. 

Appendix A 

Evaluation of Mean Free Path 

The number of collisions per unit time d© of a sample 
molecule moving at -velocity V s through target molecules in 
velocity .volume space dv5 'for hard sphere molecules is 
given (ref. 8) as 

d© = pf t V R S dvf 

where Vr is the velocity of the target molecules relative 
to the sample molecule velocity before collision and can be 
written as 

Vr = [(Vit - Vis) 2 + (V 2 t - V 2s ) 2 + (V 3 t “ V 3s ) 2 ] 1//2 (A2) 

If we subtract the local mean velocity u-j_ from both 
the velocity of the sample molecules and the target mole- 
cules we can .then write the collision rate for a Maxwellian 
distribution of target molecules as 




r B ay it <3v 2t av 3t 


(A3) 


Then transforming -variables by a rotation using Eularian 
angles (refs. 8 and 10) yields 



(A4a) 


(A4b) 


(A5a) 



13 


RAREFIED GAS DYNAMICS 
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(A5b) 




where Y r = (Vj_ 2 + V|) 1/ and ' V’ = (Y 2 + V|) 1/ Then, 
equation (A3) "becomes 

d0 M = 5/2 5 ( eXp " T") V lt" + ( V 2 t- V s) 2 + V 3t 1/2 <^t 3 (A6a) 

mjt ' C \ C j L. — 

Transforming into spherical coordinates 

V-it“ V t cos 0 sin cp, Vg t = V|. cos- cp, "^t^t sln 6 sin 9 

(A6b) 

yields 

d®M = C g ^2 [gxpC-p 2 )]^ sin cp dcp d0 dp t (A7) 

mrc ' 

where q is the nondimensionalized velocity V’/ c an(3 - 
where 

= (pi 2 + n 2 - 2q t (i s cos cp) 1 / 2 (A8) 

To obtain the total collision rate for the target mole- 
cules over all velocities, equation- (A7) must "be integrated 
over cp, 0, and- Integrating over 0 from 0 to k 

gives 


= -^^[exp(-M|)[]n|n R ^A d Mt 


^R, A 


+ M|/3q s if q g > ^ 
M- t + M-g/3d t if |J- t > 


(A10) 


Integrating equation (A9) over p-i. from 0 to °o, the 
total collision frequency 0^ is obtained for a sample 
molecule moving at velocity q s through a Maxwellian gas : 



(All) 


% “ 



Sarefied gas nmAMics 

exp(-ni) . , ^ 

' ( er %s)(^s 



The nondimensional mean free path is then found, as 
discussed in ref. 6, by dividing the sample molecule veloc- 
ity by the total collision rate, " which gives 


% V s V2 KhV s 

% C ^{? 7 ^[ eXp( " w2s) ]' * (er ^s)(^- ^)j 


(A12) 


If there were two different Maxwellian gases in the 
increment % + and f M „, the total collision rate would be 
given by © = + %£-> -where 0]y[+ and ©jyi_ are given as 

in equation (A12) but with different densities and mean 
thermal velocities. Since A s ’/D is given by Vg/©, 


. A 


s 


D 


%+ %_ 
V s ‘ V s 


1 



1 

A m _ ■ 


(A13) 


Appendix B 

The Target-Molecule Collision Partner 

The velocity distribution of target -molecule collision 
partners for a sample molecule passing through target mole- 
cules in a Maxwellian distribution (ref, 6) is 


d0 M J 

m *» 

exp(-|i|) ^ 

W M 

sin (p dcp dty. 


jexp(-p|) 

jr 3 / 2 (4rfp s ')| 

(" 1. \ 
^ 2P S ) 


The distribution in 9 for the target -molecule collision 
partners is 


Ve * f 

2it 


or 6 = 2rtR 


9 


(B2) 


The distribution of |Xj. for the target -molecule colli- 
sion partners is obtained from the marginal distribution 
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f a 


4 

jexp.(-p|) 

(a t |i RA dtJ 't ■ 

|exp(«(j|)j 

+ jr.l/ 2 (erf p g ) | 

r , i\ 


(B3) 


Then ^ can he picked from this distribution, using 

r(^ s ,d t ) + e(^ s ,u t ) 


v = 


(BA) 


1 

2 ^s , 


jexp(-|i^)J + it 1 / 2 (erf p g ) ^|i g + 
where, when p s > p^., e = 0 and 

3 

r(p s ,Pt) = -2p s M t [exp(--n2)J+ p g Vk(erf j^) - | “[ ex P(-^)] 

s 

- ~Jexp(-!i 2 )j + (erf p^.) (B5) 

and, when \± s < p t , rCdg^) = T(d s ^P s ) and 


e = -2 


jexp(p^)J^. + ld-p^ + ^ 2 exp(-|i 2 )J + ij (B6) 


To find cp for the target molecule, the distribution of 
target molecules can be written as a product of a marginal 
■times a conditional distribution 


d0(p t ,cp) 

0— = f(n t )f(T|D t ) 


(B7) 


Then it is necessary to pick cp from the conditional dis- 
tribution: 


f(9lD + ) = 


p R sin cp dcp 


2p 


(B8) 


R,A 


Picking cp for a given 


is found as before from 



16 


RAREFIED gas dynamics 


R m = 
9 


»R - K • M sl 3 


(A + |i 2 + 2|i U ) 3 / 2 " [ u_ - (i 
+■ « *b s 1 + <= 


(B9) 


After picking cp, p^., and 0, it is necessary to find the 
velocity components V^., using equation (A6b). Then, from 
equation (A4), 


T,* = V" U V” fes) - V'' fAsAs 1 + u. 


It 


rs , 


^t y s 


3t\ V V' 
v rs s 


v,--v« felv (^) rF-^' 
v 2t - htl— ) + v 2tl v; j - T 3t( TO,n 


rs 


rs s 


V 3t = 0 + ^2tf^J + V 3t(jf ' 


(BiOa) 

(BlOh) 

(BlOc) 


APPENDIX C 

Collisionless Solution 

For negligible molecular collisions, the local distri- 
bution in the channel is given by 


P f oo = 




(Cl) 


V° 


The local density is then given by p(l) = p = p + g + P m1 * 
Since there is no net flow across the channel, one can write 

P <V = "I^[ p +,0 C w,0 " ‘ p -l C w,l] °* Combining gives 



CO 1 + t 


- 1/2 > 
~w,l 



C2) 


The flow parallel to the channel wall is given by 
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U 1 “ < V 1> = u i o + — u i i 

1 p p- l,w,l 

The local temperature is given by 


(C3) 


rrT’ = o (V? + V 2 + V 2) _ i jyj2 w hi c h w hen integrated 

■w,0 3C 1 2 3 1 


w,o 


gives 


t = M 2 4-flM 2 - yA+-^t(l+ 1 1 / 2 ') (C4) 

3\ P v,0 p v,l p ^ vr,lj v ; 

The shear stress given by P = can he integrated 

'to give 


t 1 / 2 ] 


12 P+ 


^\\o " P (1 V ' Vi 1 


(C5) 


The total energy transferred given by = £{v ,2 V 2 ) + u i p i2 


is 


W 2 c 


l T 


p < 2 E gV )5/2 


fl +,0 P -1 , V 2 . p +,0 4.0 Pl -1 M w,l 

~ " t -1 + o “ , A 


p XI p 
APPENDIX D 


s St 1 / 2 

V, 1 


(C6) 


Continuum Solutions 

.laxwell's equatioh of transfer, which describes the rates 
of change of a quantity Q in a volume element of space, is 
given by 

^-(p(V i9 »=|^ (Bp 

where /iQ denotes the rate change in some molecular quantity 
Q due to molecular collisions. This term is zero where Q 
is conserved during collision. This is the case for Q equal 
to 1,V, or V 2 because of conservation of mass, momentum, 
and energy, respectively: 

- — (pu ± ) = 0 (D2) 

ox^ 
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pui 




- »< T i Y p = - 4 p i3 


Sx i 


sr ( ^ )+ £***«) ^ MY i v P )+ ir^ u’- 0 


(D3) 

(D4) 


Introducing the gas conductivity, viscosity and P as the 
static pressure gives the well known Ravier-Stokes equations: 



I p< v ; y ' 2 > - * - - « J* 


Su 


+ 2 a- - 3 6 ij 


S.u f , 


(D5) 

(D6) 


Since, in the present model, there is no x-^ or x 3 
dependence, these results simplify to give u^ = g(x 2 ), 
ft = constant across the channel^ and = -p.(5u^x ) 

(a constant), and also q 2 + u^ P 12 = ^(a constant). 


For hard-sphere molecules, assuming a small perturba- 
tion from equilibrium of the local distribution, the con- 
ductivity and viscosity can be written (ref. 6) as 


K = K^t 1 / 2 


h = 



■k ' - 75r rn_ R g C 0 
K 0 2 172 

640 (2 *) ' 

5mC~ 

!a 0 = - 

1 


(D7) 

(D8) 


With both walls stationary, q = -^(ST/Sx^) . This can be 
readily solved to give t^_Q = £l - ^1 - 

Since the pressure across the channel is constant using 
the ideal gas law, 


1/2 

( p \ _ 1 + ~kw,i + ^Wjl 

v P a Jm= 0 3(t)7 


; M=0 


(D9) 


The heat transfer across the channel can be found equa! to 





(D1C 
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When the upper wall is moving, the energy equation 
can he written as 


Un 


P 


( C p rr + 2 j K ^ 


(Dll) 


This can he rewritten using (pC^//c) = 2/3j = (5/2)R^ as 



, 4_ 2 

t + i 5 M - 15 P c 
12 0 


(D12) 


When the upper wall is moving at a velocity u - ^ and hoth 
walls are at the same temperature, integration can he per- 
formed to obtain 

t = _L M(M t - M) + 1 
15 w,l 

The equation /for : - Ppg - can he rewritten as 


= -u.tVs 


P 12 ~^0 


(D13) 

(D14) 


This can he integrated to give 




+ 15 


-1 


5^1 


- 2M 


V< ls " +M w,i> 1/2 
i 1S )l/2 Vl, . ( <l +15) " 


- sin 


- 1 / 


Mw,! 


Vi 


sm 






Vp***,!^ 


\( 15+ 4 ,i) 1/2 y 



(D15) 


This gives an implicit form of the velocity 'profile 
across the channel. This result can he used in conjunction 
with equation (D13) to obtain the temperature distribution 
across the channel. 

Similarly integrating equation (D14) from x = 0 to 1, 
the following can he obtained: 

-n 1/|/2 Pic> (u I’m 2 _+iF^r / M 


pc; 


12 _ 5]tKn/ \.i ( ' IyI fr,l +15 ) 

2 + 2 -/te 


0 


sm 


-1 


w,l 


V(M,, i+15) 


l72> 


(D16) 


To obtain the total heat transfer equation (Dll) is integrated: 
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JtA _ / ^ 1//2p 12\ M w ? l 

P V \ pC o J 2 

Appendix E 
Confidence Limits of Results 



(D17) 


For some number of sample molecules, a resulting quan- 
tity Q is obtained. . Repeating 'this process K times 
will give a sample set of Q. The average of the K val- 
ues would give Q, which is the average value of the. sample 
set. If an infinite number of values of Q were obtained, 
this would give a population of Q with some average value 
(Q) . The problem is finding the value of 6 such that the 
absolute difference between (Q) and ^ is less than 5 
with, a 95 percent probability. The central limit theorem 
states that the distribution of sample means f(Q) approaches 
a normal distribution with a mean value equal to the popu- 
lation mean and with a variance ■ equal to the population vari 
ance divided by K, the size of the sample averaged to ob- 
tain the sample mean. Since the 95 percent - confidence inter- 
val Is defined as twice the standard deviation, 6 = 2S/(/c)-*-/ 2 , 
where S 2 is the approximation to the variance of the popu- 
lation obtained by finding the- variance of the sample set 


S 2 = 


*? - 2 
(Q 2 - Q )K 

K - 1 


(El) 


This gives the 95 percent confidence interval 8 as 

n 

r (Q) - Q|<8 = 2P — — 

\K - 1 


(E2) 
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Figure 2. - Temperature distributions; wall-temperature ratio, 0.25; wall Mach number, 0. 
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Figure 4. - Dimensionless heat flux; wall-temperature ratio, 0.25; wall Mach number, 0. 
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Figure 5, - Dimensionless heat flux; wall-temperature ratio, 0.25; wall Mach number, 0 
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Figure 9; - Dimensionless total energy transfer and shear stress; wall-temperature ratio, 1; 
wall Mach number, 4. 




